Effects of growth conditions on optical quality and surface morphology of InGaAsBi
Li Jia-Kai1, 3, Ai Li-Kun1, Qi Ming1, Xu An-Hui1, Wang Shu-Min1, 2, †
State Key Laboratory of Functional Materials for Informatics, Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, Shanghai 200050, China
University of Chinese Academy of Sciences, Beijing 100049, China
Department of Microtechnology and Nanoscience, Chalmers University of Technology, Gothenburg SE-41296, Sweden

 

† Corresponding author. E-mail: shumin@mail.sim.ac.cn

Abstract

The effects of Bi flux and pressure of AsH3 on Bi incorporation, surface morphology and optical properties of InGaAsBi grown by gas source molecular beam epitaxy are studied. It is found that using relatively low pressure of AsH3 and high Bi flux can strengthen the effect on the incorporation of Bi and increase its content linearly with Bi flux until it nearly reaches a saturation value. The result from Rutherford backscattering spectroscopy (RBS) confirms that the Bi incorporation can increase up to 1.13%. By adjusting Bi and As flux, we could improve the surface morphology of InGaAsBi sample. Room temperature photoluminescence shows strong and broad light emission at energy levels much smaller than the InGaAs bandgap.

PACS: ;81.05.Ea;
1. Introduction

The host III–V compound material which a small quantity of bismuth (Bi) atoms is incorporated into has attracted much attention due to its unique properties of the strong enhancement of spin-orbit splitting and large band-gap reduction.[15] Bismuth is the largest and heaviest group V element with its isoelectronic energy level that resides in the valence bands of most III–V materials,[6] therefore, substituting As and/or P with a small amount of Bi into InP, GaAs, or InGaAs will result in anomalously narrow bandgaps[7,8] due to valence band anticrossing (VBAC).[912] Engineering the band structure for potential optoelectronic and electronic device applications is possible because of these interesting properties.[1316] Also, InGaAs alloy, which is widely used as a fundamental material in photonic devices, is expected to be turned into a promising temperature-insensitive semiconductor InxGa1−xAs1−yBiy by Bi incorporation.

In 2005, the first InGaAsBi quaternary alloy was grown on InGaAs buffer layer on an InP substrate by Gan et al. using molecular beam epitaxy (MBE).[12] Good crystalline quality was achieved with 2.5% Bi, which was confirmed by x-ray diffraction (XRD). Two years later, the same research group investigated temperature dependence of Bi behavior during MBE growth.[11] They achieved a maximum of 6% Bi concentration at a growth temperature as low as 260 °C. The growth at a high temperature of 450 °C, normally for InGaAs epitaxy growth, revealed no detectable Bi incorporation but the decrease of surface roughness from 1.067 nm to 0.328 nm.[17] In 2017, Zhou et al. demonstrated that up to 7.5% of Bi incorporation was achieved in InGaAsBi,[18] which was the highest Bi content until now. For the applications of InGaAsBi in electronic and optoelectronic devices, a short-wave infrared detector with quaternary InGaAsBi as an absorption layer and the bismuth content of about 3.2% has red-shifted the 50% cut-off wavelength from about to at room temperature.[19]

In this work we systematically study the effects of growth condition, including Bi flux and pressure of AsH3 on Bi incorporation, surface morphology and optical properties of InGaAsBi grown by the gas source molecular beam epitaxy(GSMBE) method.

2. Experimental methods

All InGaAsBi samples were grown on (100) semi-insulating InP substrates by a V90 GSMBE system. Elements In, Ga, Bi, and As2 coming from arsine (AsH3) at 1000 °C were used as source materials. Elements In, Ga, and Bi fluxes during growth were controlled by adjusting the respective effusion cell temperatures. Prior to growth, by slowly ramping up the manipulator temperature, the substrate surface was deoxidized under P2 pressure until the reflection high-energy electron diffraction (RHEED) pattern showed a clear and abrupt transformation to a (2 × 4) surface reconstruction. This surface reconstruction was adopted as a means for calibrating the substrate temperature. After that, InGaAsBi epilayers with a thickness of about were grown at a temperature as low as 300 °C measured by thermocouple for efficient Bi incorporation. Bi cell temperatures of 0, 510, 520, 530, and 540 °C were adopted. To study the influence of the As flux on the incorporation of Bi into InGaAsBi, the pressure of AsH3 was changed from 270 Torr to 360 Torr (1 Torr = 1.33322 × 102 Pa) with the Bi cell temperature unchanged. For simplicity, the InGaAsBi samples with varying Bi flux were renamed as a1, a2, a3, a4, and a5 as shown in Table 1. The InGaAsBi samples with varying pressures of AsH3 were renamed as b1, b2, b3, b4, and b5 as shown in Table 2.

Table 1.

Material growth parameters and measured properties of the grown InGaAsBi at different Bi effusion cell temperatures.

.
Table 2.

Material growth parameters and measured properties of grown InGaAsBi at different pressures of AsH3.

.

After growth, the Bi compositions were determined by Rutherford backscattering spectroscopy (RBS) with 2.275-MeV 4He2+ ions. The structures of all samples were characterized by ω–2θ scans along the (004) direction, performed in a Philips X’pert high resolution x-ray diffractometer (HRXRD) equipped with a Ge (220) four-bounce monochromator using Cu–Kα1 (λ =0.15406 nm) radiation. Their surface morphologies and roughness values were measured by atomic force microscopy (AFM). Photoluminescence (PL) and absorption measurements were conducted in a Fourier transform infrared (FTIR) spectroscopy by using a liquid nitrogen-cooled InSb detector and a CaF2 beam splitter. A diode-pumped solid-state (DPSS) laser (λ =532 nm) was used as an excitation source for PL measurement.

3. Results and discussion

The RBS is a powerful method of determining the concentration of heavy atoms in a crystal lattice composed of light elements. A typical result of sample a2 is shown in Fig. 1. The peak/step-like signal of the element Bi is observed at channel number around 480. This signal is strong evidence for Bi incorporation. From the simulations, the Bi concentration is calculated to be 0.7% for this sample. The high resolution x-ray diffraction (HRXRD) ω–2θ scans of samples a1–a5 are shown in Fig. 2. The narrow diffraction peaks correspond to InP (004) diffractions, and the relatively broad peaks correspond to InGaAsBi epilayers. It is clear that the mismatch among diffraction peaks of InGaAsBi and InP increases with the Bi content going up, which indicates that the lattice constant increases with the addition of Bi. As shown in Table 1, increasing Bi content brings about a larger full width at half maximum (FWHM) due to the increase of Bi-related dislocations in the alloy of sample.[20] From the peak separation between InP and InGaAs1−xBix in the XRD spectrum, each lattice constant can be estimated. For convenience, the variations in In composition and lattice deformation have been ignored.[21] Based on Vegard’s law, the lattice constant (C) of the quaternary alloy In0.53Ga0.47As1−xBix can be expressed as follows:

Fig. 1. (color online) RBS spectrum from the InGaAsBi film with a Bi concentration of 0.7%. The black, violet, red, blue, magenta, olive, and navy lines represent the randomly and simulated spectra of the InGaAsBi films with the individual P, Ga, As, In, and Bi, respectively.
Fig. 2. (color online) HRXRD (004) ω–2θ scans for InGaAsBi epilayers on InP with Bi content of 0% (a1), 0.7% (a2), 0.98% (a3), 0.99% (a4), and 1.13% (a5).

Lattice constants are assumed to be 6.639, 6.234, 6.058, and 5.653 Å for InBi, GaBi, InAs, and GaAs in zinc blende structures, respectively.[22] Therefore, the amounts of Bi in the InGaAsBi alloys estimated from the lattice constant are around 0%, 0.70%, 0.98%, 0.99%, and 1.13%. The Ga has an effect to reduce the lattice constant after being incorporated into InP and InAs lattice.[23] It is expected that the lattice constant of InGaAsBi can be matched to that of InP by appropriately controlling the Ga and Bi content.

The surface morphologies are characterized by AFM scan using tapping mode for samples a2–a5. The root-mean-square (RMS) roughness versus the Bi content of the InGaAsBi epilayer is shown in Table 1. Figure 3 shows AFM images of the samples. All the samples contain some holes. The formation of these holes is likely to be due to thermal evaporation of oxides from the InP substrates. While the Bi concentration increases to 1.13% in sample a5, big Bi related metallic dots can be seen in the fourth image of Fig. 3, which is in the center of white rectangle. We suggest that Bi surface migration plays an important role. The rough starting surface of InP substrate after oxide desorption restricts the Bi surface diffusion, resulting in a local increase of Bi concentration. So a smooth surface facilitates Bi surface diffusion and thus uniform Bi incorporation.[24]

Fig. 3. (color online) Surface morphologies of the InGaAsBi samples a2 (a), a3 (b), a4 (c), and a5 (d).

Figure 4(a) shows the PL spectra of the samples a1–a5 at room temperature. The InGaAs reference sample (sample a1) shows a weak PL signal at 0.74 eV which is close to the bandgap of In0.53Ga0.47As[25] and originates from the band-to-band (BB) transition. After a small amount of Bi is incorporated (sample a2), a strong and broad PL spectrum appears between 0.78 eV and 0.66 eV, peaked at about 0.715 eV. As expected, the PL peak energy decreases with increasing Bi content, implying that the band gap of InGaAs1−xBix alloydecreases with increasing Bi content. The peak energies are 0.714, 0.691, 0.688, and 0.682 eV for the samples with 0.70% (a2), 0.98% (a3), 0.99% (a4), and 1.13% (a5) Bi, respectively. A bandgap reduction coefficient of 48 meV/Bi% is deduced and calculated from the experimental data. The PL intensity decreases dramatically when the Bi content increases from 0.7% to 0.98%, which implies that non-recombination centers in In0.53Ga0.47As1−xBix increase because of the deterioration of crystalline quality. The FWHMs of the PL spectra of the samples a2–a5 are nearly 0.060, 0.064, 0.065, 0.075 eV, respectively, which confirms that the deterioration in crystalline perfection of InGaAsBi is due to the increasing Bi content. Figure 4(b) shows the variation of the FWHM of the PL spectra of samples a2–a5 with Bi content for samples a2–a5. The trend implies that the deterioration process in crystalline perfection of InGaAsBi with higher Bi content becomes quicker than with relatively low Bi content. It is because at higher Bi effusion cell temperature, more Bi atoms will be accumulated on surface, forming large metallic Bi droplets on the growth front surface, which hinders the InGaAsBi from uniformly growing.

Fig. 4. (color online) Room-temperature PL spectra for the samples with different Bi content. (b) Variations of the FWHM of the PL spectra of samples a2–a5 with their Bi content.

The InGaAsBi samples, with pressures of AsH3 varied but other conditions fixed, are renamed as b1, b2, b3, b4, and b5 as shown in Table 2. The HRXRD ω–2θ scans of samples b1–b5 are shown in Fig. 5. Through the simulations and comparison of the HRXRD ω–2θ scans of samples b1–b5 with the HRXRD ω–2θ scan of sample a2, the amounts of Bi in samples b1–b5 are estimated from the lattice constant at about 1.05%, 0.48%, 0.21%, 0.10%, and 0%. It is obivious that the lattice constant increases with the pressure of AsH3 decreasing. As is well known, both As and Bi are group-V atoms which tend to bond with group-III atoms. So an excess As supply will outcompete Bi since III–As bonding is stronger than III-Bi bonding. Hence, the decreasing of AsH3 pressure will spontaneously enhance the possibility of III–Bi bonding and thus increase the Bi content. Then, the diffraction peak in Fig. 5 shifts toward the small-angle side with the pressure of AsH3 decreasing.

Fig. 5. (color online) HRXRD (004) ω–2θ scans for InGaAsBi epilayers (b1–b5) on InP at different AsH3 pressures.

Figure 6(a)6(e) show AFM images of b1–b5 samples. These images reveal that as the pressure of AsH3 decreases from 360 Torr to 300 Torr, the root mean square (RMS) roughness values of samples decrease from 0.361 nm to 0.196 nm. After the pressure of AsH3 reaches 300 Torr, the RMS roughness value begins to increase, which peaks at 3.210 nm when the pressure of AsH3 is 240 Torr. Figure 6(f) shows the variation of the RMS roughness values of samples with the pressure of AsH3. This trend and the relationship between RMS roughness value and pressure of AsH3 can be qualitatively comprehended like this: both As and Bi are group-V atoms which tend to bond with group-III atoms. Hence, the decreasing of AsH3 pressure will increase the Bi content. So when the AsH3 pressure decreases from 360 Torr to 300 Torr, this trend will facilitate the incorporation of Bi atoms so that fewer excess Bi atoms will stay or accumulate on the surfaces of the samples. Then, the RMS roughness values of samples decrease. However, the continuous decreasing of AsH3 pressure will not result in a monotonic increase of Bi content, but leads to forming droplets because excess group-III atoms on the growing InGaAsBi surface cannot be evaporated. So when the AsH3 pressure continously decreases from 300 Torr to 270 Torr, the RMS roughness values of samples increase from 0.196 nm to 0.311 nm because the excess group-III atoms on the growing InGaAsBi surface cannot be evaporated at such a low growth temperature. The RMS roughness value of sample peaks at 3.210 nm when the pressure of AsH3 is 240 Torr.

Fig. 6. (color online) Room temperature PL spectra for the samples b2–b5 at different AsH3 pressures.

The PL spectra at room temperature of b2–b5 samples at different AsH3 pressures are shown in Fig. 7. As the AsH3 pressure decreases, the PL peak is red-shifted, which confirms that the lower AsH3 pressure enhances the incoporation of Bi atoms. The PL intensities show their small differences at high AsH3 pressures ranging from 360 Torr to 300 Torr but decrease dramatically at about 270 Torr AsH3 pressure and we cannot observe the PL peak of sample b1, which implies the deterioration of the crystalline quality because the excess group-III atoms on the growing InGaAsBi surface cannot be evaporated from the surface.

Fig. 7. (color online) (a)–(e) Surface morphologies of InGaAsBi samples b1, b2, b3, b4, b5. (f) Variations of the RMS roughness values of the InGaAsBi samples b1, b2, b3, b4, b5 with AsH3 pressure.
4. Conclusions

In this work, we investigate the effects of Bi flux and AsH3 pressure during growth on Bi incorporation, surface morphology and optical properties of InGaAsBi grown by GSMBE on InP substrate. It is found that using high Bi flux and low AsH3 pressure during growth can strengthen the effect on the incorporation of Bi. The Bi content increases linearly with Bi flux until it reaches nearly a saturation value. The maximum incorporated Bi content is found to be 1.13%. A relatively low Bi flux and a moderate AsH3 pressure during growth can improve the surface morphology. We obtain the photoluminescence spectra of InGaAsBi at room temperature, which show strong and broad light emission at near infrared wavelength. The low AsH3 pressure will remarkably damage the optical properties of InGaAsBi grown by gas source molecular beam epitaxy.

Reference
[1] Oe K 2002 J. Appl. Phys. Jpn. 41 2801
[2] Kudrawiec R Kopaczek J Misiewicz J Petropoulos J P Zhong Y Zide J M O 2011 Appl. Phys. Lett. 99 251906
[3] Lu X Beaton D A Lewis R B Tiedje T Zhang Y 2009 Appl. Phys. Lett. 95 041903
[4] Fluegel B Francoeur S Mascarenhas A Tixier S Young E C Tiedje T 2006 Phys. Rev. Lett. 97 067205
[5] Carrier P Wei S H 2004 Phys. Rev. 70 035212
[6] Alberi K Dubon O D Walukiewicz W Yu K M Bertulis K Krotkus A 2007 Appl. Phys. Lett. 91 051909
[7] Wang K Gu Y Zhou H F et al. 2014 Sci. Rep. 4 5449
[8] Gu Y Wang K Zhou H F et al. 2014 Nanoscale Res. Lett. 9 24
[9] Lewis R B Masnadi-Shirazi M Tiedje T 2012 Appl. Phys. Lett. 101 082112
[10] Sweeney S J Jin S R 2013 J. Appl. Phys. 113 043110
[11] Feng G Oe K Yoshimoto M 2007 J. Crystal Growth 301 121
[12] Gan F Masahiro Y Kunishige O Akiyoshi C Yuji H 2005 J. Appl. Phys. 44 L1161
[13] Jin S Sweeney S J 2013 J. Appl. Phys. 114 213103
[14] Marko I P Batool Z Hild K et al. 2012 Appl. Phys. Lett. 101 221108
[15] Dongmo P Zhong Y Attia P et al. 2012 J. Appl. Phys. 112 093710
[16] Gu Y Zhang Y G Song Y X et al. 2013 Chin. Phys. 22 037802
[17] Wang L Zhang L Yue L et al. 2017 Crystals 7 63
[18] Zhou S Qi M Ai L Wang S Xu A Guo Q 2017 J. Appl. Phys. Jpn. 56 035505
[19] Gu Y Zhang Y G Chen X Y Ma Y J Xi S P Du B 2017 Appl. Phys. Lett. 108 032102
[20] Zhao H Malko A Lai Z H 2015 J. Crystal Growth 425 89
[21] Petropoulos J P Zhong Y Zide J M O 2011 Appl. Phys. Lett. 99 031110
[22] Oe K Asai H 1996 IEICE Trans. Electron. E79C 1751
[23] Wang K Wang P Pan W W et al. 2015 Semicond. Sci. Technol. 30 094014
[24] Gu Y Zhang Y G Li A Z Wang K Li C Li Y Y 2009 Chin. Phys. Lett. 26 077808
[25] Vurgaftman I Meyer J R Ram-Mohan L R 2001 J. Appl. Phys. 89 5815